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Abstract 



■ Spectrum of density perturbations in the Universe generated from quantum- 

gravitational fluctuations in slow-roll-over inflationary scenarios with the 
Brans-Dicke gravity is calculated. It is shown that after inflation the isocur- 
vature mode of perturbations may be neglected as compared to the adiabatic 
mode, and that an amplitude of the latter mode is not significantly different 
from that in the Einstein gravity. However, the account of the isocurvature 



, mode is necessary to obtain the quantitatively correct spectrum of adiabatic 

perturbations. 

in 

o> 

1 I. INTRODUCTION 

cr 

Inflationary expansion in the early Universe not only is able to explain the observed 
degree of homogeneity, isotropy and flatness of the present-day Universe but can also 
account for the origin of small initial density perturbations which produce gravitationally 
bound objects (galaxies, quasars, etc.) and the large-scale structure of the Universe [fj. 
Historically, among models of inflation making use of a scalar field (called an inflaton), 
the original model or the first-order phase transition model M failed due to the graceful 
exit problem, which was taken over by the new || and the chaotic [Q inflation scenarios 
where the inflaton scalar field is slow rolling during the whole de Sitter (inflationary) stage. 
The latter property was shared by the alternative scenario with higher-derivative quantum 
gravity corrections (where the role of an inflaton is played by the Ricci scalar R) just 
from the beginning. Note that a simplified version of this scenario - the R + R 2 model 
- was even shown to be mathematically equivalent to some specific version of the chaotic 
scenario J7| (see also a review in ||). In order to obtain a small enough amplitude of density 
perturbations in all these slow-roll-over models, the inflaton should be extremely weakly 
coupled to other fields. It is therefore not easy to find sound motivations to have a such a 
scalar field in particle physics. (See, however, [[| for recent improvements along this line.) 



*To be published in the proceedings for the fourth workshop on general relativity and gravitation 
edited by K. Maeda, T. Nakamura, and K. Nakao 



1 



Reflecting such a situation, the extended inflation |H| scenario was proposed several 
years ago to revive a GUT Higgs field as the inflaton by adopting non-Einstein gravity 
theories. Although the first version of the extended inflation model, which considers a 
first-order phase transition in the Brans-Dicke theory |11|| , resulted in failure again due to 



the graceful-exit problem ||12|| , it triggered further study of more generic class of inflation 



models in non-Einstein theories, in particular extended chaotic inflation |13| where both the 
inflaton and the Brans-Dicke scalar fields are in the slow rolling regime during inflation. 
Note that the natural source of Brans-Dicke-like theories of gravity is the low-energy limit 



of the superstring theory [n,r5j with the Brans-Dicke scalar being the dilaton. 

Several analyses have been done on the density perturbations produced in extended new 
or chaotic inflation models fl6H20ll, all of which made use of the constancy of the Bardeen's 



gauge-invariant quantity ( or its equivalent on super-horizon scales and match it directly 
to quantum field fluctuations at the moment of horizon crossing which would be the correct 
procedure in a single component inflationary model. However, such a treatment is not correct 
in the present case which contains two sources of quantum fluctuations, namely, the inflaton 
and the Brans-Dicke scalar field. In such a model not only adiabatic but also isocurvature 
perturbations are generated and ( does not remain constant for all modes. 

In the present article we analyze density perturbations generated in both new and chaotic 
inflationary models in the Brans-Dicke gravity. Keeping the above point in mind and ex- 
tending the method used in |22|,|23| to find spectra of all modes of adiabatic and isocurvature 
fluctuations in the multiple inflationary scenario in the Einstein gravity (with scalar fields 
interacting with each other through gravity only), we separate adiabatic and isocurvature 
modes at the inflationary stage and carefully evaluate their final amplitude. As a result, 
we arrive at a formula which is apparently very different from that used previously but find 
that constraints on the model parameters in both models remain approximately the same 
as in Einstein gravity if we express them in terms of the values of coupling constants at the 
end of inflation, contrary to the analysis by Berkin and Maeda . 

We make use of the conformal transformation which transforms the original or the Jordan 
frame to the Einstein frame in which equations are somewhat simpler. We calculate the 
spectrum of fluctuations in the Einstein frame and then interpret the result to the Jordan 
frame. 

The rest of the present article is organized as follows. In §2 we define the Lagrangian and 
perform conformal transformation to write down the equations of motion in the Einstein 
frame. Then in §3 solutions to these equations are given and they are interpreted in the 
Jordan frame. They are applied to both chaotic and new inflationary models in §4. Finally 
55 is devoted to discussion and conclusion. 



II. BASIC EQUATIONS 

A. Conformal Transformation 

It has been shown that with the help of a conformal transformation a wide class of 
non-Einstein gravity models can be recast in the action 
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S = J C^-gd'x = J LLfc + \g^d, X d,X - U(x) + ^^d^d^ - e^W^) 



(1) 



where k 2 = 8irG, (3 and 7 are constants, \ and are the dilaton and the inflaton fields, 
respectively ||16|| . For example, in the case of the Brans-Dicke theory |TT[, the original action 
in the Jordan frame 



S 



+ —^g^d^ BD d^ BD + V<9,<R0 - v{<t>) 

167T 167T<I>BD 2 



#6? X, 



(2) 



is transformed into (|l|) through the conformal transformation 



— O 2 " 



n 2 



K 



57T 



BD 



exp 



fu + 3/2 



(3) 



with /3 = 27 = , 2 and = 0. Observations constrain u to be u > 500 pi], so that 
y/uj+3/2 — 

(3 < 0.09 <C 1. Since the Brans-Dicke field $bd remains practically constant in the radiation 

or matter dominant stage (see below), x must be equal to zero at the end of inflation in 

order to reproduce the correct value of the gravitational constant today. Although we focus 

on the original Brans-Dicke theory in the present paper, we will not put the explicit values 

for (3 and 7 below until the very end of the calculation in order to sustain applicability of 

our analysis to other non-Einstein gravity theories. 



B. Background equations in the Einstein frame 

Since the analysis is much simpler in the Einstein frame ([I]), we first investigate equations 
of motion derived from ([!]) and then transform the results into the original frame. Taking 
the background as the spatially flat Friedmann-Robertson- Walker spacetime, 

ds 2 = dt 2 - a(t) 2 dx 2 , (4) 

the scalar field equations for the homogeneous parts read 

X + 3H X + - PKe-^W^) = 0, (5) 

4> + 3H<j) - 7 kx</> + e^-^W (0) = 0, (6) 

with an overdot denoting time derivation. Because (3, 7 <C 1 and V(<p), being the inflaton's 
potential, realizes slow roll over of 0, the above equations can be approximated as 

3Hx = PKe-P K W(<l>), (7) 
3#0 = -e^ )KX V'{(f>), (8) 



with 



d\ 2 k 2 



tf 2 =(r) =y^vm (9) 
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during inflation. We note that the above approximation is valid provided that the following 
inequalities are satisfied. 

max{e- 7KX 2 ,x 2 } < e^ KX V{(/)), 

|V'(0)| < 3/€e- 7Kx/ V(0), (10) 
V"{<f>)e^® KX < H 2 . 



Then the system admits a solution as a function of the scale factor as fl"3| . [T6 



X=-lna = --z, (11) 

K K 

2 r<t> V(cp) , l-e~^ z 

where a subscript / denotes the value of each quantity at the end of inflation and we have 
set df = 1. 

In the case that <fi changes slowly compared with x an d that V(0) can be regarded as a 
constant, evolution of the scale factor and % is given by 

a(t)=a t^, x(t) = ^lnf j. (13) 

We thus have a power-low inflation with an extremely large power exponent. 



C. Linear perturbations 

We now turn to linear perturbation, which is most easily taken into account in a gauge- 
invariant manner using the longitudinal gauge, 

ds 2 = (1 + 2§)dt 2 - a(i) 2 (l - 2*)6 ij dx i dx> , (14) 

where $ and ^ are gauge-invariant variables related with Bardeen's $^ and as $ = $^4 
and \I/ = —&h p5 |. Assuming an e* kx spatial dependence, each Fourier mode satisfies the 



following equations of motion which are derived from the perturbed Einstein equations. 

$ = ^, (15) 

$ + = y (x^X + e"^<Mx), (16) 



+ 3#5x + "' 



* _ l2^L e - 7KX 2 + (/3«) V^W) ] 5x + i^ KX 

a 2 2 1 



-i3Ke- 0KX V'{(j))8(j) = 2{x + 3#x)$ + $X + ^X + 7Ke" 7KX 2 $, (17) 



4 



+ (3H - 7«*)<ty + e (l -P )KX V"((f)) ) 50 - 7 /t05x 



a 2 



+ (7 - [3)KV'((f))e^-' J)KX 5x = 2(0 + 3#0)$ + $0 + - 2 7 k;0x$, (18) 

where 8\ and 50 are gauge invariant fluctuation variables of the respective fields and we have 
suppressed the argument k. Under the slow-roll approximations fllCf ) the last two equations 
read 

S'x + 3H5 X + ( ^ + (^) 2 e~^ KX V((j)) ) 5 X + iKe^ KX 4)8^ - j3Ke~^ KX V' 



= 2/^e- 7KX \/(0)$ + $x + 3^x, (19) 

50 + 3#50 + + e^-^W'i^ 50 - 7 /t#x + (7 - ^V'^e^^x 

= -2e (7 - /3)KX 1/'(0)$ + $0 + 3^x- (20) 

Since what we need are the non-decreasing adiabatic and isocurvature modes on large 

scale k <C aH, which turn out to be weakly time-dependent as will be seen in the final result 

p3| , we may consistently neglect $ and those terms containing two time derivatives. Then 
the equations of motion are greatly simplified to 

* = + e"^050) = ^S X - ^50, (21) 

3H8'x + (Pn) 2 e- pKX V(<P)5x - ^e'^W (<P)5<P = 2/^ e - /3KX K(0)$, (22) 

3H6'(f> + e^ K *V"(0)50 + (7 - f3)KV\<P)e^-P )KX 5 X = -2e { ^ )KX V'{(j))^. (23) 

III. GENERATION OF PERTURBATIONS 
A. Solution of non-decreasing modes 

We shall now solve eqs. (|2"ID-(p3|) . Inserting ( pi]) to (|22"D we easily find 3H5x = 0, or 



Similarly (|23|) is transformed to 



S x = P-Q-l = const. (24) 



3H5(j> + e^~® KX (V(0) - Y^-j H + /57e (7 " /3)KX V r/ (0)Q 1 = 0. (25) 

The above equation can be solved taking as an independent variable instead of t. With 
the help of (§) we find 
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3 4^4=- eM, " xr <' (26) 



so that fl2"5] ) reads 

Its solution is given by 



6<f> + p-yQi. (27) 



u r + ^ = - w, (28) 



^ 1/(0) V UfV'{if) ^ K 2 J K 2 V((f)) 

with Q2 being another integration constant. We therefore find 

Thus we have obtained a generic solution of the system (l2l"|)-(|23"l) containing two undeter- 
mined constants. 

In order to clarify physical meaning of the above solution, we should divide it to adiabatic 
and isocurvature modes. The latter mode is characterized by its vanishingly small contri- 
bution to the gravitational potential $, while the growing adiabatic mode can be described 
by the following universal expression. 

= £ /V)*^, (31) 

x 4> a Jo 

where C\ is a constant and the latter approximate equality in each expression is satisfied 
during the inflationary stage. 

From ®-(D we find ' 

x = ^ * = ~ H ^vW' lP = Y + ^{vm) ' (32) 

and it turns out that defining new constants C\ and C3 by Qi = C\ — C3 and Q2 = — C3 
discriminates between adiabatic and isocurvature modes in the final result. That is, 

(33) 
(34) 

(35) 





G\ C3 




X 


H H' 




5(f) 

J~ 


§ + §(e"^-l), 


$ = 







V v{4>) J 2 



In the above expressions, terms in proportion to C\ and C3 represent adiabatic and isocur- 
vature modes, respectively. The isocurvature nature of the C3— terms is guaranteed by the 
fact that the second term in the right-hand-side of (|35|) is vanishingly small in the last stage 
of inflation when we have x — 0. 
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B. Quantum fluctuations 



We shall next determine the constants C\ and C3 from amplitudes of quantum fluctua- 
tions of the scalar fields generated during the inflationary stage. Thanks to the inequalities 
(|10|), eqs. (|19|) and ( p0|) can be approximated by equation of motion of a free massless 
scalar field in inflating background for k > aH and even in the region k < aH but with 
H(tk) 3> \H(tk)\(t — tk) where tk is the time k— mode leaves the Hubble horizon during infla- 
tion. The standard quantization gives the well-known result, that is, the Fourier components 
of the fields can be represented in the form 



;(k), ty(k) 



H(t 



'2k 3 



e^(k), 



(36) 



where the exponential factor in the latter equality is present because <fi has a non-canonical 
kinetic term in the action in the Einstein frame. Here e x (k) and e^(k) are classical random 
Gaussian quantities with the following averages. 



(e x (k)> = fo(k)) = 0, (e,(k)e*(k')> = 5^(k- k'), i,j = X A- 
We thus find 



c 3 



e< KX H 5 -i + (1 - e 7KX )#— 
X 



H\t k ) 



V X, 



H\t k ) 
V^k 3 



2k 3 

rt(k) - 



g2 7 K X I _ e 7«X 

' — e^(k) + e x (k) 



X 



' - ' 1 • ' x ( k ) 



X 



(37) 

(38) 
(39) 



Note the interesting fact that in the specific case of the "soft" inflation [gg] (7 = 0), the 
Brans-Dicke dilaton does not contribute to the adiabatic mode (Ci) at all. 



C. Adiabatic modes 



From ( |3~0"D curvature perturbation due to primordially adiabatic fluctuation is given by 

-1 



ad 



1 + 



3(l + w) 



n - P 

Cx, w = -, 

V 



(40) 



with p and p being total pressure and energy density, respectively, since $ is related to the 
gauge-invariant density fluctuation Sp^'/p, which reduces to 5p/p in the comoving gauge, 
as 



c) 



2 /_fcj 

3 \#a 



we find 



5p 



(<0 



./ 



e^H 5 -t + (1 - e 7KX )#— 
X. 



(41) 



(42) 



7 



at the horizon crossing (k = aH), where / is a constant equal to 4/9 during radiation 
domination and to 2/5 during matter domination (we omit the minus sign here because it 
may be absorbed into the stochastic variables 8(f> and 8\)- 

The above formula is quite different from that used by the previous authors 116-2011, 
who did not account for the isocurvature mode and its mixing with the adiabatic mode 
during inflation, and concluded that a density fluctuation at the second horizon crossing is 
approximately equal to 8p/(p + p) evaluated at the first horizon crossing during inflation. 
That is, up to a factor of order of unity, 

8p H(\x\8 X + e~^\4>\8<f)) 
— ~ : (incorrect), 

P X 2 + e~^(f> 2 



(43) 



which gives 



5 -l ~ { f > l0|e " /2 (incorrect) (44) 



P 



Note also that we would have reached the above formula ([43]) or its analogue in the Jordan 
frame fl2"0"| , had we put C 3 = neglecting the isocurvature mode in fl35f) . Clearly this is not 
justifiable because it results in imposing a redundant constraint between 8\ and 8(f) in eqs. 
([33]) and ([341). The difference between the expressions (§^) and (f4~3"D is especially large near 
the points where V'((f>) = 0. Then Eq. ( f42"|) predicts much larger adiabatic perturbations 
than Eq. §|). 

We are interested in the amplitude of density fluctuations about 40 ~ 60 e-folds before 
the end of inflation, or at z ~ 40 ~ 60, corresponding to large-scale structures within the 



horizon scale today. From (|il|) the corresponding value of x ranges kx = —40/3 ~ —60/5, 
for which we have e 7KX ~ 1.2 with u ~ 500. Thus the above incorrect formula and ours give 
practically opposite results with respect to the two cases |x|> \4>\ and |x|< \<f>\- In some 

models of inflation, unless f3 is too small, we find |^| < so that the adiabatic fluctuation 
is given by 

5 -P * e™™*, (45) 
P 101 

which gives practically the same result as in the Einstein gravity. On the other hand, in the 
limit u — > oo or (3, 7 — > 0, we find e 7KX = 1 so that the second term in (|42|) does not 
contribute to adiabatic fluctuation and we again find 

8p H8(f) 



P ' 101 ' 

in agreement with the result in the Einstein gravity. 



(46) 



D. Isocurvature modes 

Isocurvature fluctuations can be significant if some ingredient of the model is essentially 
decoupled from the usual matter and constitutes a part of dark matter. In our model, the 
Brans-Dicke dilaton is a candidate of such an ingredient. 



S 



During inflation, its fractional comoving energy perturbation satisfies the inequality, 
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^^ = 4(^)_„ = <*(« <c , (47) 

Px Px+Px dt\ X J ^ \ V((p) J 

Since C3 is of the same order of C\, which gives the amplitude of adiabatic perturbation at 
the second horizon crossing, fractional density fluctuation in x is smaller than the adiabatic 
one. Furthermore, its contribution to the total energy density fluctuation, 8p^ / p turns out 
to be much smaller than the adiabatic counterpart, because we find p x = 0{uj~ 1 )p <C p in the 
post inflationary universe as will be seen below. We, therefore, conclude that isocurvature 
fluctuations are negligible in this model. 



E. Density fluctuations in the physical (Jordan) frame 

We can relate the amplitude of density fluctuations in the Einstein frame with that in the 
Jordan frame using the conformal transformation f27fl. First we note that gauge-invariant 
variables in both frames are related with each other as 



$ = $ 



5Q 

TP 



* = * + 



5Q 



where 5Q = ^Sx^i is gauge-invariant perturbation of Q. Using the equalities 



a = Q 1 a, di = u 1 dt, H 



one can write (|30|) and ([31]) in terms of physical variables as 

1 



1 



H + -± 



VtJ aVL 2 



atfdt' 



$X 
X 



n 



Xt 



Ci 
dQ 



dn 2 dt'. 



We, therefore, find 



Of 



VtJ dVL 2 



dtt 2 dt 



1 



H 

dtf? 



dn 2 dt' 



(48) 
(49) 

(50) 
(51) 

(52) 
(53) 



The above two equations are the desired formula to relate density fluctuations in the Einstein 
frame and those in the Jordan frame through C\. In practice, however, since the Brans- 
Dicke field varies extremely slowly in the post-inflationary universe, one can regard Q as a 



constant in these equations. For example, in the matter dominated stage, we find |TT 

' i 



3u>+4 



$BD(*) = $BD0 



to, 



a(t) = d [ j- J 



(54) 
(55) 
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Hence Q^/Q is smaller than H by a factor less than 1/500. Thus we may simply conclude 
that the adiabatic fluctuations are described by the same formula in both frames. 

Using (|54D , we can estimate energy density of x m the Einstein frame in this stage as 

1. 2 2oo + 3 1 „ 1 
P - = 2 X = W+~W^ = 6^ (56) 

Thus we may conclude the isocurvature perturbation due to x is n °t important as stated in 
the previous subsection. 



IV. APPLICATION TO SPECIFIC INFLATIONARY MODELS 

Having developed general considerations, we now apply the above results to two specific 
inflation models, namely chaotic inflation || and new inflation |4[], and examine constraints 
on their model parameters as well as the spectral shape. 



A. Chaotic inflation 

Here we consider two different potentials V((j>) = ^m 2 2 and V((f)) = which are 
denoted collectively as V((f>) = — (f> n . In this model, inflation occurs at large and it is 
terminated when \4>/(f>\ becomes as large as H at <p — <Pf — \fnj k. Since we are interested 
in the last 60 e-foldings in the inflationary stage, we may approximate (12) as 



K 



dip 



K 

2n 



n 



z. 



(57) 



The error in the last expression is only about 12% with z = 60, (3 = 2^ = 0.09. The e-folding 
number, Zk, when the comoving wave- number k leaves the Hubble radius during inflation 
satisfies 



k_ 



2z k ) 



for 1 < z k < 60. 



(58) 



We can express the amplitude of curvature perturbation on scale I = ^ as 



i + 



i+ 



3(1 + w) 
2 

3(1 + w) 



2fc 3 (|C 1 | 2 ) 



2tt 



- 1 2 
K 



2ttV 3n 



-Zk 



n 



(59) 



again for 1 z k < 60. 

Since the large-angular-scale anisotropy of background radiation due to the Sachs- Wolfe 
effect is given by ST/T = $/3, we can normalize the value of X n by the COBE observation 
28[. For = 2 7 = 0.09, we find 
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8T 1 4/ . f 9^- for n = 2 

_ = 60) s ( ^ forn = 4 (60) 

= 1.1 x 10" 5 (61) 

Note that our normalization of the background solution a/ = l, x/ = ( see Eqs. (p^JT^ )) 
means that m and A4 are values of physical parameters ("coupling constants") at the end 
of inflation. However, since \ generally grows only as logarithm of t after inflation and is 
even constant during the whole radiation-dominated stage in the Brans-Dicke case f3 = 27, 
their present values are not significantly different from those at the end of inflation. So, it 
is natural to use just the latter values when comparing non-Einstein gravity models having 
variable coupling constants with models based on the Einstein gravity. We find 

m = 1 x 10 13 GeV, n = 2, (62) 
A 4 = 1 x 1(T 13 , n = 4 (63) 

which is no different from the values obtained assuming the Einstein gravity [2~9~|. Since the 
behavior of the system approaches to that in the Einstein gravity as we increase u, we can 
conclude that in Brans-Dicke theory, model parameters of the inflaton's potential should 
take the same value as in the Einstein gravity. 



B. New inflation 



Next we consider new inflation with a potential 



V(<j>) = V - ^0 4 , (64) 



for which we find, from ([12]), 

* Viif) ^^K 2 V fl 1 \ K 2 



K 2 



vW) dip ^^-^J^^- (65) 



We can again express the amplitude of curvature fluctuation as a function of z^, which is 
now related with k as, 



l r e(^h H , (66) 



$(/) 



1 + 



3(1 + w) 



-1 



e 2 



ij8G8-7)*fc, 



!\ „ p~~ Zk 



(67) 



with H f = yjf V . Taking (3 = 2 7 = 0.09 again, it predicts the amplitude of ST/T to be 
compared with COBE data as 



ST H 
— (z k ~ 60) -21^ + 0.4— L 
1 Mpi 



(68) 
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Since Hf should also satisfy 

inr-< 10~ 5 (69) 

M P r y 1 

to suppress long-wave gravitational radiation of quantum origin |3(J , we find 

A< 2 x lCT 13 (70) 

from (|6q). Again its amplitude is practically no different from the case of the Einstein 
gravity. 



V. CONCLUSION 



In the present paper we have investigated density perturbations generated during infla- 
tion in the Brans-Dicke theory, paying attention to the fact that quantum fluctuations of 
both the inflaton and the Brans-Dicke dilaton fields contribute to them. This implies that 
not only adiabatic (curvature) but also isocurvature perturbations exist and that the simple 
conservation formula of a perturbation outside the Hubble radius cannot be used. 

Working in the Einstein frame, which is connected with the original frame through a 
conformal transformation, we have obtained expressions for growing adiabatic and isocur- 
vature modes and have shown that only the former is important after inflation. Then we 
have corrected the previously used formula based on the use of incorrect matching of super- 
horizon perturbations with quantum fluctuations at horizon (Hubble radius) crossing. Since 
both the inflaton and the dilaton turn out to contribute to the perturbation with the same 
order of magnitude, the amplitude itself is not significantly altered compared with the case 
of inflation in the Einstein gravity. 

However, our constraints on the model parameters turned out to be more severe than 



those obtained by Berkin and Maeda | lq| . This is because they have taken the initial value 



of x rather arbitrary, which might be appropriate to the case of original soft inflation model 



in which inflaton's potential is multiplied by an exponential potential by hand [26|. On the 



other hand, we, concentrating on the Brans-Dicke model, had to normalize its value to be 
zero after inflation to reproduce the correct gravitational constant today. 

In conclusion, values of inflaton's coupling parameters in the Brans-Dicke gravity (if 
taken at the end of inflation) which are required to produce the correct amplitude of present 
adiabatic perturbations are practically no different from those in the Einstein gravity in the 
case of both chaotic and new inflationary models. 
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